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mentum distribution of Higgs bosons produced via gluon-gluon fusion at 
the LHC. In particular we focus on the joint resummation formalism which 
resums both threshold and transverse momentum corrections simultane- 
ously. A comparison of results obtained in the joint and the standard 
recoil resummation frameworks is presented. 
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1. Introduction 

The Large Hadron Collider (LHC) is widely expected to be the Higgs 
boson discovery machine. In particular, the dominant channel for the pro- 
duction of a light Higgs particle at the LHC is gluon fusion gg — » HX 
The search strategies for the Higgs boson rely deeply on the knowledge of 
production characteristics, the transverse momentum (Qt) of the produced 
boson being one of the most important quantities. In this talk we describe 
an application of the joint resummation formalism |2j to calculate transverse 
momentum distribution of Higgs bosons produced through the gluon fusion 
mechanism at the LHC. 

It is a general feature of perturbative calculations in QCD that close to a 
phase space boundary partonic hard-scattering cross sections acquire large 
logarithmic corrections. These corrections are related to soft and collinear 
gluon emission and arise from cancellations between virtual and real con- 
tributions at each order in perturbation theory. The threshold and recoil 
corrections are the two notable examples often discussed in this context. 
The threshold corrections of the form a™ ln 2ri_1 (l — z)/(l — z) become large 
when the partonic cm. energy approaches the invariant mass Q of the pro- 
duced boson, z = Q 2 /s — > 1. The recoil corrections, in turn, are of the 
form a™ In 2 ™ -1 (Q 2 / 'Q T ) and grow large if the transverse momentum carried 
by the produced boson is very small, Qt <C Q- Thus, sufficiently close to 
the phase-space boundary, i.e. in the limit of soft and/or collinear radia- 
tion, fixed-order perturbation theory is bound to fail. A proper treatment 
of higher-order corrections in this limit requires resummation of logarithmic 
corrections to all orders. 

In the Standard Model the leading 0(a 2 ) process for Higgs boson pro- 
duction via gluon-gluon fusion proceeds through a heavy quark loop, with 
the top quark loop providing the most significant contribution. In the limit 
of large top mass, m t — > oo, the Higgs coupling to gluons through the top 
loop can be described by the effective ggH vertex jSj ■ This simplifying ap- 
proximation was shown to be valid up to a few percent accuracy in the case 
of NLO calculations [I]. The fixed-order predictions for the total production 
rate are currently known at the NNL (next-to-next-to- leading) order. Al- 
though not as large as the NLO corrections [3], the NNLO corrections were 
found to be substantial, increasing the NLO predictions by around 30% jH). 
Moreover, it was shown that the prevailing contribution to these corrections 
corresponds to the soft and collinear gluon emission [51 |7j , thus reinforcing 
the need for a careful treatment of logarithmic corrections to all orders. 

The resummation techniques are well established both in the thresh- 
old [SI [5] and in the recoil ^| case for the Drell-Yan production-type 
processes. The Drell-Yan mechanism and the mechanism for Higgs boson 
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production through gluon-gluon fusion are similar. This makes it possi- 
ble to apply (after implementing necessary changes accounting for gluons, 
instead of quarks, in the initial state) the already developed resummation 
methods to Higgs boson production. The resummed predictions were ob- 
tained in [IJ ^] for threshold resummation and in |13L I14L 1151 116j for recoil 
resummation. A joint, simultaneous treatment of the threshold and recoil 
corrections was first introduced in [2J • It relies on a novel refactorization 
of short-distance and long-distance physics at fixed transverse momentum 
and energy Similarly to standard threshold and recoil resummation, 
exponentiation of logarithmic corrections occurs in the impact parameter b 
space Fourier-conjugate to transverse momentum Qt space as well as 
in the Mellm-iV moment space (HI El) conjugate to z space. The resulting 
expression respects energy and transverse momentum conservation. A full 
phenomenological analysis of Z boson production at the Tevatron in the 
framework of joint resummation can be found in 18 i , whereas Higgs boson 
production at the LHC was studied in [TU] . 



2. The jointly resummed cross section 

The general expression for the jointly resummed cross section applied 
to Higgs boson production via gluon fusion, reads [T§] 
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where r = Q 2 /S, rrih is the mass of the Higgs boson, and irT<jQ5(Q 2 
denotes the lowest order partonic cross section for the process gg 
the limit of large mj, with 

h _ V2G F a 2 s (m h ) 
a °~ 576^ ' 
The function H contains the hard virtual part of the perturbative cor- 
rections and up to 0{a s ) is given by pfl l2"T] 
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Apart from a lower limit of integration, at the next-to-leading-logarithm 
(NLL) accuracy the form of the Sudakov factor E^J(N, b,Q,n) in the jointly 
resummed expression is the same as for the recoil resummation |18j : 
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The functions A and B are perturbative series in a s and their coefficients 
can be determined by comparing fixed order predictions with an expansion 

of resummed result |211 120 j . NLL accuracy requires using A ( g \ B^ l) and 

Af ] in Eq. 0: 
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where Ca = 3, Cp = 4/3, Tr = 1/2, and Np is the number of flavors. The 
higher order (needed at NNLL) coefficient B 2 is also known |21| 



(6) 



The quantity x(iV, b) appearing in the lower limit of integration in © 
is specific to joint resummation. The LL and NLL logarithmic terms in the 
threshold limit, N — > oo (at fixed b), and in the recoil limit b — > oo (at fixed 
iV) are correctly reproduced with the following choice of the form of x 



x(N,b) = b + 



where r/ is a constant and we define 



N 



1 + Tfb/N' 



(7) 



(8) 
0) 



b = bQe^/2 , 

with jp the Euler constant. 

The functions C(Q,b, N, /i, fj,p) in Eq. Q are given by: 



C a /H(Q,b,N,n,LL F ) = ^2C a/j (N,a s (n)) S jk (N,Q/x, Hf) fk/H{N,^ F ) ■ 

j,k 

(10) 

The product of parton distribution functions fk/H a ^ scale fj,p, and a matrix 
can be seen as corresponding to parton densities evaluated at the scale 
Q/X- The evolution from the scale \x F to Q/x is accurate to NLL in x and 
represented by the matrix £ (N,Q/x, ^f)- The origin and the structure 
of the evolution matrix £ was discussed in detail in Ref. |18l I19j . The 
coefficients C a /j(N,a s ) have a structure of a perturbative series in a s , and 
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are determined in the same way as for recoil resummation, i.e. up to 0(a s ), 
C 9/g (N,a s ) = 1 + g^l + g^, (11) 
C 9/q {N,a s ) = = gcV^L_ = C g/q (N,a s ) . (12) 

The expression is formally accurate to the next-to-leading-logarithm 
(NLL) level. However, due to the large colour charge of the incoming gluons, 
the cross section exhibits increased sensitivity to the Sudakov logarithms. 
It is known for recoil resummation that the NNLL terms have a significant 
impact on numerical results. Motivated by this finding, we decide to include 
the NNLL terms containing the coefficient, in the way consistent with 
the recoil resummation. In our formalism we also include the 0(a s ) per- 
turbative expansions for the functions H and C, which formally give NNLL 
contributions. As discussed in Ref. the NNLL resummed cross sections 
are resummation scheme dependent, and the choice of the resummation 
scheme is reflected in the value of the coefficients B^ 2 \ and We 
exercise the freedom of the resummation scheme choice by demanding that 
the function H, calculated with a s taken at the scale Q, collects the hard 
virtual part of the NLO corrections. The Sudakov factor and the C coeffi- 
cients contain then only soft or collinear contributions. The values of 
and C^ 1 ' listed above are for this particular resummation scheme. 

By incorporating full evolution of parton densities the cross section (^Q) 
correctly includes also the leading a™ \n 2n ~ 1 (N)/N collinear non-soft terms 
to all orders. At the NLO, this can be seen by expanding jointly resummed 
cross section to 0{a s ) (here for illustration purposes only in the gg channel), 
integrated over Qt 



& gg = a (0) ^£ | _ ACa ln 2 N + g7r6o ln at + ll + 3vr 2 
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where Si(N) = J2j=i J 1 = i>(N + 1) + 7_e, with tp the digamma function. 
In the large iV limit this gives 

*» = °f> £ Uoa m 2 N + 4C A ^L + ii + 37r 2 1 + Q feN\ | (M) 

which can be compared to the the large iV limit of the NLO result in the 
Mellin space 
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The agreement between the expanded jointly resummed expression and the 
exact NLO result down to the 0(1/N) is clear. The mismatch in the con- 
stant 7r 2 term between ()14p and ()15|) originates from the value of the Cg) g 
coefficient and is a NNLL effect. A development of the joint formalism at 
the NNLL would eliminate this disagreement as well as provide a way to 
include other NNLL coefficients, most notably B^ 2 \ 

2.1. Numerical results and discussion 

In order to have predictive power, the resummed expression Q needs 
to be supplemented by a definition of the inverse Mellin and Fourier trans- 
forms from N and b space. In the joint approach, the inverse integrals are 
both treated as contour integrals in the complex space of N and b. For the 
integrals to be well defined, the contours must not run into the Landau pole 
or singularities associated with the form of the function \. This procedure 
provides an unambigous definition of resummed perturbation theory with- 
out an introduction of additional dimensional scales and implies a functional 
form of non-perturbative corrections. We refer the reader to Ref. [T2I f° r a 
detailed discussion of the parametrization of the contours. 

The joint resummation formalism with the inverse transforms defined 
as contour integrals ensures that predictions can be obtained for any non- 
zero value of Qt- This is not possible in the standard recoil approach 
without adding some non-perturbative term of the form —gb 2 to the Sudakov 
exponent. However, for the purely technical reasons of numerical stability 
we also include such a factor in our jointly resummed cross section. The 
value of the g parameter, g = 1.67 GeV 2 , is adopted from the study in 
Ref. [21]. However, we checked that the dependence of the results at small 
Qt on the value of g is negligible, in agreement with what was found for the 
case of pure b space resummation. At large Qt, where the ln(Q T / Q 2 ) terms 
taken into account by resummation lose their importance, it is necessary 
to match the resummed result with a fixed-order result. Here the jointly 
resummed result is matched to the 0(a s ) perturbative result ^3], in the 
way described in [T8] . 

The numerical results for the Higgs boson transverse momentum dis- 
tribution calculated in the joint resummation framework were obtained as- 
suming nih = 125 GeV, (i = fj,p = Q = and using CTEQ5M [23] parton 
distribution functions. The parameter n in the definition of \ J7J) is chosen 
to be n = 1. We checked that the numerical dependence of the predictions 
on the value of n is small |19j . 

Apart from the joint resummation predictions, in Fig. 1 we also show 
the recoil-only (i.e. \ = b) resummed result. At small to moderate Qt, the 
b space resummed prediction is slightly higher and broader than the one 
provided by the joint resummation but the difference is small. Consequently, 
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we conclude that the threshold effects are of modest importance at these 
values of Qt and the pure recoil resummation is fully applicable there. 

The Q-r-integrated joint cross section, by definition, is expected to return 
the threshold resummed result. Although it does so formally up to NLL, 
numerically the integrated joint distribution returns a result which is ~ 10% 
lower than the threshold cross section. We find that this suppresion is caused 
by subleading terms included in the joint resummation, more specifically 
terms oc l/(iV — 1) in the expansion of the joint expression (j!3|) . These 
terms arise from our treatment of evolution in the coefficients C, cf. Eq. (|1U|) . 
They are important only in the small N limit and therefore not present in 
the threshold resummed expression. (However, it is interesting to observe 
that the NLO cross section in Mellin space contains the same subset of 
terms cx 1/(N — 1) as the expanded joint expression taken at b = 0, plus an 
additional subset of terms oc l/(A r — l) 2 . The numerical effects of the two 
subsets cancel almost entirely, leading to a relatively good approximation 
of the NLO cross section by the threshold resummed result.) The small N 
limit corresponds to the limit of small z = Q 2 /s <C 1. Given partonic cm. 
energies available at the LHC, the small z terms of the form a™ ln^ 2n_1 ^ (z) / z 
may indeed play a significant role for light Higgs production. These terms 
can be resummed on their own |25| . but their full inclusion in the joint 
formalism alongside the threshold and recoil corrections requires further 
work. 



sqrt(s)=14 TeV, M n =125 GeV, CTEQ5M 




L i i I i i i I i i i I i i i I i ""rT^^^T I i 

20 40 60 80 100 120 

Q T (GeV) 

Fig. 1. Transverse momentum distribution for Higgs production at the LHC in the 
framework of joint resummation and of "pure-QT" resummation. 
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